Steady-state MreB helices inside bacteria: dynamics without motors 
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Within individual bacteria, we combine force-dependent polymerization dynamics of individual 
MreB protofilaments with an elastic model of protofilament bundles buckled into helical configu- 
rations. We use variational techniques and stochastic simulations to relate the pitch of the MreB 
helix, the total abundance of MreB, and the number of protofilaments. By comparing our simula- 
tions with mean-field calculations, we find that stochastic fluctuations are significant. We examine 
the quasi-static evolution of the helical pitch with cell growth, as well as timescales of helix turnover 
and denovo establishment. We find that while the body of a polarized MreB helix treadmills towards 
its slow-growing end, the fast-growing tips of laterally associated protofilaments move towards the 
opposite fast-growing end of the MreB helix. This offers a possible mechanism for targeted polar 
localization without cytoplasmic motor proteins. 

PACS numbers: 87.16.-b,87.16.Ac,87.16.Nn 



I. INTRODUCTION 

The eukaryotic cytoskeleton organizes cell shape, cell 
polarity, cell division, and non-diffusive subcellular trans- 
port. F-actin, microtubules, and intermediate filaments 
comprise the cytoskeleton, and act with the associated 
proteins that provides spatial and dynamic control of cy- 
toskclctal function [1]. Prokaryotic cells have cytoskeletal 
analogues, such as the FtsZ-ring associated with cell di- 
vision [2] together with its "divisome" of associated pro- 
teins. Bacteria also have a number of polymerizing cy- 
toplasmic proteins, such as ParM [3] and MinD [2] that 
exhibit distinctive helical structures within the cell. 

Recently, the actin homologue MreB has been shown 
to play a cytoskeletal role in many bacteria [2, 4, 5]. 
MreB forms a continuous cytoplasmic helix that runs the 
length of nearly all rod-shaped prokaryotes, including Es- 
cherichia coli, Bacillus subtilis and Caulobacter crescen- 
tus [6] , and it has been implicated in shape determination 
and polar protein localization. 

In most Gram-positive bacteria MreB is present to- 
gether with several paralogues, such as Mbl and MreBH 
in B. subtilis. The helical pitches for MreB or Mbl, sep- 
arately observed by immunoflourescence microscopy, are 
reported to be 0.73 ± 0.12 /im and 1.7 ± 0.28 /im, respec- 
tively [7]. More recent measurements of fluorescent fu- 
sions of MreB and of Mbl report pitches of 0.6 ± 0.14 /im, 
with colocalization [8]. In Gram- negative species, such 
as E. coli and C. crescentus, only MreB is present. In 
E. coli, pitches of 0.46 ± 0.08 ^m have been reported [9]. 
In all cases, the helices are dynamic, with elements mov- 
ing along the main helix at reported speeds ranging from 
6nm/s [10] to 70nm/s [11]. The helical structure has 
also been observed to condense into a ring at midcell 
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near the time of division in E. coli [12], C. crescentus 
[13] and B. subtilis (where only MreBH coils) [14]. 

The MreB helix appears to be composed of a bun- 
dle of individual "protofilaments" [15-17]. Quantitative 
immunoblotting has been used to measure the molecu- 
lar abundance of various MrcBs. In B. subtilis, there 
are roughly 8000 MreB monomers and 12000 - 14000 
monomers of Mbl [7] while E. coli has roughly 17000 — 
40000 monomers of MreB [9] . Neglecting the cytoplasmic 
fraction of monomeric MreB, these abundances suggest 
a bundle thickness of about 10 protofilaments [18]. 

In B. subtilis Mbl is necessary for proper insertion of 
new peptidoglycan, which occurs in a helical fashion [6], 
while MreBH is necessary for the localization and func- 
tion of the cell wall hydrolase LytE that is believed to 
recycle the outer layers of the cell wall, also in a helical 
fashion [14]. Cells with mutant mreB are wide, rounded 
and usually not viable [19]. Helical bundles of MreB may 
contribute to the spatial localization of associated MreC, 
MreD, and PBP2 that, in turn, help to coordinate cell 
wall synthesis. It has also been suggested that helical 
filaments of MreB paralogues under tension can lead to 
spiral morphologies [20]. 

Disruption of MreB leads to loss of proper polar lo- 
calization of a number of proteins such as the chemo- 
taxis protein Tar and the Shigella flexneri virulence fac- 
tor IcsA in E. coli [21], and three integral membrane 
proteins (PleC, DivJ, CckA) in C. crescentus [22]. Polar 
localization in C. crescentus was disrupted by either un- 
derexpression or overexpression of MreB. When normal 
MreB expression was returned, polar localization was re- 
established [22] . This suggests that MreB has a continual 
role in either direct polar trafficking of these proteins or in 
the maintenance of landmarks necessary for their proper 
positioning [4]. 

The polar proteins in C. crescentus are normally di- 
rected towards distinct (stalked and swarmer) poles in 
different stages (swarmer, stalked, and predivisional) of 
its life cycle. For example, PleC is localized to swarmer 
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poles in swarmer and predivisional cells, DivJ is local- 
ized to stalked poles, while CckA is localized to both 
poles of predivisional cells [23]. After MreB expression 
is disrupted and restored, PleC and DivJ are restored 
randomly to either pole [23]. This suggests that MreB 
may be polarized within C. crescentus and that the po- 
larity of the MreB helix is randomly restored after its 
disruption. However, tracking of individual YFP-MreB 
molecules shows unpolarized motion [10], raising ques- 
tions about the mechanism of specific polar targeting. 
MreB-directed targeting to specific poles has not been 
reported in E. coli or B. subtilis. 

MreB interacts with both RNA polymerase (RNAP) 
[23] and SetB, a chromosome defect suppressor [24], and 
has been implicated in the fast polar translocation of the 
origin-proximal regions (oriC) [27] of newly-rcplicatcd 
DNA in C. crescentus [25] and in E. coli [23] (see however 
[26]). Time- lapse microscopy has shown that the polar 
transport of oriC in B. subtilis had an average speed of 
2.8nm/s and a peak speed of 4.5nm/s [27]. 

MreB is a homologue of the eukaryotic cytoskeletal 
protein actin [28, 29] . Actin filaments are used to change 
cell shape and to move bacteria such as Listeria mono- 
cytogenes via polymerization forces [30], and in muscle 
contraction and for organelle movement via collections 
of myosin motor proteins [1]. Actin assembly is regu- 
lated through a number of "actin-binding proteins" that 
variously control cross-linking, bundling, filament nucle- 
ation, end-capping, filament cutting, monomer seques- 
tration and desequestration. MreB, in contrast, does not 
have any clearly identified motor proteins or associated 
proteins that regulate polymerization. Notably, MreB fil- 
aments spontaneously bundle in vitro without associated 
proteins [31]. 

The varied roles of MreB inside bacterial cells raise 
some basic questions. What is the origin of its helical 
configuration, and how significant are the forces that the 
MreB helix applies? What does the helical pitch of MreB 
filaments depend on? What aspects of the MreB system 
can be understood in terms of actin-like polymerization 
dynamics? Specifically, must we invoke yet-to-be dis- 
covered prokaryotic motor proteins or accessory proteins 
controlling MreB polymerization to explain MrcB-relatcd 
polar localization of proteins such as Tar, IcsA, DivJ, 
PleC, and CckA? Finally, does the small size of the bac- 
terial cell affect MreB polymerization, as compared to 
actin polymerization in much larger eukaryotic cells? 

To address these questions, we present a model of the 
MreB helix with stochastic polymerization dynamics of 
protofilaments together with global elasticity of a heli- 
cal MreB bundle constrained by the bacterial cell wall. 
Our model provides a quantitative relationship between 
helical pitch, total abundance of MreB protein in a par- 
ticular cell and the thickness of the protofilamcnt bun- 
dles. The bundled MreB protofilaments are in a dynam- 
ical steady-state, and undergo constant advection as the 
polymerized subunits treadmill. We discuss how this ad- 
vection could be harnessed for targeted polar localization 



without motor proteins. The steady-state dynamics also 
allow us to address other dynamical processes such as 
protein turnover in FRAP experiments [32] or recovery 
from A22-induced disruptions of the MreB helix [22] . 



II. MODEL 
A. Protofilament polymerization 

Both actin and MreB polymerize into polarized fila- 
ments. Addition and dissociation of monomers occur 
at the ends of the asymmetric filament, both at the 
"barbed" ("+" or fast-growing) tip and the "pointed" 
("— " or slow-growing) tip. The kinetics of actin 
polymerization are well-characterized by concentration- 
dependent polymerization rates k^ n c and fc£ n c at the 
barbed and pointed ends, respectively, where c is the 
cytoplasmic monomer concentration, and concentration- 
independent depolymerization rates k^ e and fc^ ff . 

If a force F is applied to a filament's tip, the polymer- 
ization rate is reduced. When thermal bending fluctua- 
tions are much faster than the polymerization dynamics, 
the polymerization rates at either end of the filament are 
reduced to [30] 
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where oq — 5.1 nm is the MreB monomer length [33] 
and fcsT = 4.1 pNnm at room temperature. This force- 
dependent polymerization rate can also be obtained from 
thermodynamic arguments in the high- force limit [34]. 
We apply it to MreB polymerization dynamics within 
the cell. 

In the absence of force, each filament grows above, and 
shrinks below, the critical cytoplasmic concentration 
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The asymmetry between the polymerization and depoly- 
merization rates at barbed and pointed ends leads to 
treaclmilling, in which the filament length remains con- 
stant while depolymerization from the pointed end is bal- 
anced by polymerization at the barbed end [1]. The force- 
independent treadmilling rate is [35] 
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If the filament position is fixed, treadmilling results in a 
net advection of all polymerized monomers towards the 
pointed end. 

The kinetic rate constants for MreB polymerization are 
yet to be determined explicitly, but appear to differ sig- 
nificantly from eukaryotic actin. In addition to sponta- 
neously nucleating and bundling with much greater ease 
than actin, purified MreB from Thermotoga maritima 



3 



polymerizes in vitro faster and exhibits a much lower crit- 
ical concentration (c c MreB = 0.003 ^<M [31] compared 
with c cact j n = 0.17 [10]). Nevertheless, in vivo ob- 
servations of single molecule motion of MreB in C. cres- 
centus suggest that the treadmilling rate is similar to 
actin [10]. By starting with in vitro rate constants from 
eukaryotic actin [29] and scaling the on-rates by a factor 
°f c c,actin/ c c MreB = 55, we recover the reported critical 
concentration of in vitro MreB from T. maritima. Fol- 
lowing in vitro observations that the treadmilling rate 
is AM re B = 1.2s -1 [10], we further scale all four rates 
by ^MreBMactin = 2.1 — preserving the MreB tread- 
milling rate. These scalings preserve the pointed/barbed- 
end asymmetries of actin and represent the least-intrusive 
modification of actin polymerization dynamics to make it 
consistent with observed MreB dynamics. The resulting 
barbed-end polymerization rate constant, k on (see Ta- 
ble 1), is close to the diffusion limit [36] indicating that 
the pointed/barbed-end asymmetries of MreB may differ 
significantly from actin and/or that the appropriate c c 
for in vivo measurements may be significantly above the 
T. maritima value. However, our qualitative results do 
not depend on the precise parameter values used in this 
paper. 

We model polymerization dynamics by the 
stochastic addition/dissociation of monomers at the 
barbed/pointed ends of individual MreB protofilaments 
[17] using the scaled kinetic rate constants discussed 
above and listed, together with other parameters, in 
Table 1. Except for growing cells in Sec. Ill B, a standard 
cell length of L c = 3/um and cell radius of R c = 400 nm 
are used. 



B. Bundle ultrastructure 

The ultrastructure of the MreB helix - the precise ar- 
rangement, orientation and length distribution of the in- 
dividual protofilaments that make up the helical bundle 
- remains a mystery. Several hypotheses have been put 
forward [16] and Fig. 1 illustrates five basic possibilities. 
Several of these are less plausible. The slippery arrays in 
Fig. IE are unlikely to be able to support the forces the 
cytoskelcton must withstand as it pushes against the cell 
wall, and recent biochemical experiments have demon- 
strated large lateral interactions between filaments [31]. 
The ultrastructure in Fig. 1C leads to tapered bundles as 
antiparallel protofilaments, which cannot slide past each 
other, treadmill in opposite directions. Such tapering is 
not seen experimentally. In this paper we therefore con- 
sider ultrastructures composed of polarized bundle(s) of 
protofilaments: either one bundle (Fig. 1A) or two an- 
tiparallel bundles that freely slide with respect to each 
other (Fig. ID). Since the clastic and polymerization 
properties of the second case follow straight-forwardly 
from the former, we will mostly focus on a single polar- 
ized non-slipping filament bundle (Fig. 1A) and revisit 
the possibility of antiparallel bundles slipping with re- 



spect to each other in the final discussion. 

It is quite possible that individual MreB protofilaments 
do not continuously extend from one end of the bacterial 
cell to the other, similar to actin cables in yeast [37]. In- 
deed, in C. crescentus individual protofilaments appear 
to be much shorter than the cell length, only 392 ± 23 nm 
on average [10]. For the purposes of our model, the me- 
chanical and end-polymerization properties of discontin- 
uous bundles of protofilaments are equivalent to bundles 
of continuous protofilaments. Systematic heterogeneities 
in the MreB helix thickness have not been reported, but 
our model does not depend on how the cell regulates the 
average number of protofilaments in a cross-section of the 
filament bundle. Protofilament association, dissociation 
and nucleation are thus implicitly included in our model. 
Polymerization and depolymerization away from the fila- 
ment edges may occur and do not affect our steady-state 
results: whatever the distribution of lateral exchange, 
in steady state monomer incorporation and separation 
are balanced. Of course, significant numbers of active 
protofilament tips away from the cell poles can affect the 
dynamical timescales that we discuss. 

C. Elastic bundle 

In vitro, MreB typically polymerizes into straight fil- 
aments [2]. However, MreB adopts ring-like coiled con- 
figurations in spherical mutants of normally rod-shaped 
organisms [38], and forms helices in rod-shaped cells [15]. 
These observations suggest that normally straight elas- 
tic MreB protofilaments may simply buckle into helices 
inside the cylindrical confinement provided by the rela- 
tively hard cell wall. 

A self-consistent model for the observed MreB helices 
consists of a particular ultrastructure of protofilaments, 
as in Fig. 1, buckled into a helical configuration by the 
cylindrical cell wall. A steady state exists where the poly- 
merization force at the tips of the helices is balanced 
by the mechanical forces of the helical configuration. 
We treat the cell as a spherocylinder with total length 
L C + 2R C and radius R c and total volume ttR^,(L c + 2R c ). 
The helix is assumed to extend throughout the cylindri- 
cal part of the cell, but not into the hemispherical poles, 
as indicated by experiment [7]. 

Although helical equilibria of clastic filaments have 
been investigated since the 1800's, they remain a con- 
temporary topic [39]. In the elastic Cosserat model [40], 
a filament is parametrized by its unstretched arclength 
s e [0,Lfu] where Lfu is its total unstretched length. 
The Hamiltonian of such a filament is 

n = \j^'\{B K {sf + Cr{sf)n{sf 

+ E(1- \dr/ds\fn{s)]ds (4) 

where f(s) is the position of the centcrline, k(s) is the 
local curvature, t(s) is the local twist and n(s) is the 
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local filament thickness (measured in number of protofil- 
amcnts). B, C and E are the bending, twisting and 
stretching modulii of an individual protofilament, respec- 
tively. For an actin bundle, the thickness dependence of 
bending and twisting ranges from linear (n) for slippery 
protofilaments to quadratic (n 2 ) for non-sliding protofil- 
aments, with a crosslinker dependent crossover [41]. For 
MreB, lateral interactions appear strong so we assume a 
quadratic dependence on n. 

By imposing the observed helical configuration we can 
use variational techniques to estimate the forces working 
against monomer addition at the filament tips. In Ap- 
pendix A, we derive the force acting at the filament tips 
along the filament direction in the inextensiblc (E — > oo) 
and freely twisting (r(s) = 0) regime: 



/ B sin 2 0(l + 3cos 2 0)(n 2 ) L fil > L c 
L c > Lfu, 



(5) 



where 9 is the pitch angle of the helix and fs = B /2R 2 w 
0.031 pN is the elastic force scale. L c and R c are the 
length and radius of the cylindrical portion of the cell, 
respectively. If the bundle thickness n(s) exhibits signif- 
icant inhomogeneity, then the appropriate average thick- 
ness n is the root-mean-square average thickness along 
the bundle length. Other than the buckling point at 
Lfu = L c , Fb is independent of Lfu for a given 8. The 
helical pitch p and pitch angle are related by 
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2irR c 
tan 9 



cos 9 



Lfu 



(6) 



As the pitch p vanishes, 9 — > n/2, and Fb — » jsn 2 ■ The 
force Fb has a maximum of Fg ax = ^fsn 2 at 9* = 
arctan\/2 w 0.96 w 55°. 

Many of our results reflect the fact that we work in 
a regime where elastic forces, Fb, and the polymeriza- 
tion force scale, jo = ksT/ao — 0.80 pN, are similar in 
magnitude. 



III. RESULTS 

A. Stochastic steady-state 

In a typical B. subtilis cell of volume 1.9 /jm 3 , only 
4 unpolymcrizcd monomers arc necessary to achieve the 
critical concentration of MreB, c c = 0.003 /iM — suggest- 
ing that stochastic effects, due to the discrete molecular 
nature of the polymerizing monomers, may be significant. 
Similarly, the number of protofilaments in a typical cross- 
section of the MreB helix is small (less than 20), and the 
number of protofilaments in contact with the cell wall 
at the helix tips, n t i P , is even smaller. A deterministic 
mean- field analysis of the steady-state (see Appendix B), 
neglecting stochastic fluctuations, can be compared with 
fully stochastic simulations to explore the impact of var- 
ious stochastic effects in this system. 



We simulated n protofilaments that grew and shrank 
stochastically within a common pool of No monomers 
according to the force-dependent polymerization rates, 
Eq. 1, and force- independent depolymerization rates. 
The forces on each protofilament tip were determined by 
the constraint force on the bundle, Eq. 5, divided among 
the leading filaments at that tip, where 1 < nu p < n. 
The filament bundle was studied once it reached a steady- 
state. The length of each of the n protofilaments and the 
monomer concentration continued to fluctuate within the 
steady-state, as did the number of protofilament tips, 
nup, at a given end of the filament bundle. 

For a given cell geometry (L c and R c ) and total number 
of monomers (No), each bulk thickness n yields a unique 
helical steady-state configuration with a particular pitch 
p and average cytoplasmic monomer concentration c. As 
seen in Fig. 2, as the abundance -/V increases for a given 
number of protofilaments n, the pitch decreases due to 
longer bundles. Conversely, at a given No, thicker bun- 
dles (larger n) leads to larger pitches. 

The rectangles in Fig. 2 represent independent ex- 
perimental measurements of MreB abundance and cable 
pitch in B. subtilis [7, 14] and E. coli [9]. In B. subtilis, if 
the three MreB isoforms (MreB, Mbl, and MreBH) bun- 
dle together into a triplex structure, the total number 
of monomers should be the sum from each homologue, 
which we estimate is 23000 - 40000. Wc sec that with 
current experimental pitch and abundance estimates, E. 
coli has n € [3,16] protofilaments in each polarized MreB 
helical bundle. Under the assumption that they are bun- 
dled independently of the other isoforms in B. subtilis 
we estimate that Mbl has n <G [9,21], while MreB has 
n G [3,7], with pitches from [7]. Under the assumption 
that the isoforms are mutually non-slipping, we estimate 
n G [8, 21] in the triplex with the mutual pitch from [8]. 
It is therefore possible that the distinct pitches reported 
by [7] and [8] are the result of changes in the bundling 
properties of the variously labeled and tagged MreB ana- 
logues. 

The concentration of MreB in monomeric form shown 
in Fig. 3 is low, representing less than 0.1% of the total 
cellular MreB. This is in contrast to the same model using 
eukaryotic actin kinetics, for which 1 — 5% is monomeric 
(data not shown), and in dramatic contrast with obser- 
vations of FtsZ polymerization, for which 30% is associ- 
ated with the Z-ring and 70% is diffuse in the cytoplasm 
in vivo [42]. Notwithstanding, the average monomeric 
concentration of MreB is significantly above the critical 
concentration due to a reduced polymerization rate which 
arises from the constraining force at the bundle tips. The 
shape of the curves in Fig. 3 follows from the force vs. 
pitch relationship in Eq. 5, which grows quickly, reaches 
a maximum then decreases again as the pitch angle ap- 
proaches 7r/2. As we shall see, force and cytoplasmic 
concentration increase together, while, as we have seen, 
pitch and No are inversely related. As a result, for each 
n the cytoplasmic concentration exhibits a similar max- 
imum vs. Nq as force does vs. pitch. This maximum, 
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corresponding to 9* , is at larger Nq for larger n. We also 
see in Fig. 3 that the fluctuations in the monomeric con- 
centration are very large and approximately independent 
of n and N . The cell is often instantaneously below the 
critical concentration despite the upward bias due to the 
constraint forces. 

As seen in Fig. 4a, the maximal force F sustainable 
by the bundle in steady state increases with the cyto- 
plasmic concentration c. For this part of the figure c 
was held fixed, to facilitate comparison with both the 
mean-field calculations in Appendix B and the stochas- 
tic n = 2 calculations in Appendix C. The increase of 
F with c comes from Eq. 1, the force- induced reduc- 
tion of the polymerization-rate, and reflects the excess of 
the monomer fraction over the critical concentration line 
seen in Fig. 3. F increases with the number of protofila- 
ments, n, at a fixed cytoplasmic concentration, since it is 
distributed over the leading nu p protofilamcnts. The ex- 
cess of F over the maximal mean-field prediction (dashed 
lines), where nu p — n, is recovered analytically for n = 2 
(solid line) in Appendix C. The excess arises because any 
of the nup protofilaments may grow and extend the bun- 
dle length, after which the remaining near-tip protofil- 
aments will quickly catch up without force retardation. 
The result is enhanced growth at a given F, or, equiva- 
lently, a larger F at which steady-state is reached. The 
stochastically-enhanced force is more significant as n in- 
creases, and also increases with c due to the increased 
polymerization rate. 

In Fig. 4b, we examine the additional effects on F due 
to the large fluctuations in c that were apparent in Fig. 3. 
Here we plot Fb vs. the average c for fully fluctuating 
filament bundles. Various c were explored by varying Nq. 
We include the mean-field results (dashed lines) for ref- 
erence. Two differences with Fig. 4a are apparent. The 
first is that the fully fluctuating results have a maximal 
force sustainable by the elastic bundle, Fg^ = 4/3/sn 2 . 
With c fixed we could use an arbitrarily stiff bundle to ex- 
plore a wider range of forces, but here we need to choose 
a particular stiffness (by equating F = Fb ) to accurately 
couple fluctuations in c with those of F. The second dif- 
ference is that F for the fully-fluctuating system is below 
the mean-field results for n < 5, and is lower than the re- 
sults in Fig. 4a for a fluctuating n t i p . This systematic de- 
crease can be seen as arising from averaging the concave- 
down curves from Fig. 4 over the very large c fluctuations, 
which leads to a stronger decrease for smaller average c. 

There are significant effects on the forces sustained 
by the filament bundle both due to fluctuations in the 
small number of protofilaments supporting the force at 
the bundle tip, n tip , and due to the large fluctuations in 
the monomer concentration, c. Because of the curvature 
of the F vs c curve, these fluctuation effects modify the 
mean-field steady-state bundle force in opposite direc- 
tions. While for small numbers of protofilaments (n < 5) 
the overall effect is to reduce F for a given c, the net effect 
for the physiological range of bundle thicknesses (Fig. 2) 
and cytoplasmic concentrations (Fig. 3) is an increased 
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B. Cell growth 

The MreB helix grows as the cell grows, doubling its 
lateral length before dividing. Our quasi-static approach 
can accommodate cell growth. We assume that the to- 
tal number of monomers Nq is proportional to the cell 
length and that the number of protofilaments n is length- 
independent. We simulated cell growth in a regime to- 
wards the low-iVo end of Fig. 2, with N = 3 x 10 3 L C so 
that an average 3 /zm cell contains 9000 MreB monomers, 
and towards the high-7V end, with N = 10 4 L C . Fig. 5 
shows how the steady-state helical pitch varies as the cell 
length, L c , ranges between 2 — 4/im. 

For most bundle thicknesses, n, the pitch is nearly con- 
stant as the cell elongates. However, for thicker bundles 
in the low-iVo regime, the helical pitch increases signif- 
icantly as the cell doubles in size. This increase is due 
to fluctuations at small cell lengths, L c . While it ap- 
pears strange that at large n, the stochastic effects are 
larger, the pitch is determined by the filament length 
L fu , which only depends on the maximum protofilamcnt 
length. While the mean protofilamcnt length fluctuates 
less with increasing n, the maximum protofilamcnt length 
is an extremal property of the bundle — and increases 
with increasing n. At small L c , the relative fluctuations 
in the cytoplasmic fraction also increases. For the larger 
iVo regime of Fig. 5b, the pitches are much smaller, Lfu 
larger, and relative length fluctuations correspondingly 
smaller. 

In the physiological regimes shown by boxes in Fig. 2, 
pitch should not change significantly during cell growth 
if the overall concentration of MreB monomers remains 
constant. Mbl in B. subtilis may have longer pitches in 
longer cells [11, 15], though this is not a strong effect 
[7]. Experimental observations of considerable variabil- 
ity of the number of helical turns per cell within cells 
of the same strain, size and growth conditions [11] may 
imply corresponding variability of MreB expression or of 
the cross-sectional number of protofilaments n — which 
complicates analysis. 



C. Macromolecule trafficking 

A vital role of MreB is the polar localization of proteins 
such as Tar in E. coli [21], the cell polarity markers DivJ 
and PleC in C. crescentus [22], and the origin-proximal 
regions of the newly-replicated chromosome in E. coli, 
B. subtilis and C. crescentus [23]. One possibility is 
that these passengers associate with yet-to-be-discovered 
motor proteins that use MreB as a track to the poles 
[21]. A second possibility is that these proteins simply 
bind to the helix and advect with the continuous tread- 
milling, eventually ending up at one of the polar tips. A 
third possibility is that they associate with leading tips of 
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MreB protofilamcnts [16], perhaps via intermediary pro- 
teins analogous to formin for actin filaments. Here we 
quantitatively analyze these possibilities, which are de- 
picted in Fig. 6, together with associated translocation 
speeds with respect to the fixed bacterial axis. 

For any transport mechanism, a characteristic speed 
along the filament bundle, v tr ans, translates into a speed 
v z = cos Ovtrans relative to the cell's axis. Assuming the 
protein initially binds at a uniformly random location 
along the bundle, the average time to reach a pole is 
then 

{ttrans} ^ 7; = ~Z 7. • (7) 

2v z 2cos 0v trans 

This time can be compared to the cell division time to 
see whether it provides a plausible mechanism for polar 
localization. 

The myosins that transport organelles along actin 
tracks in eukaryotes travel at speed v mot ~ 200 — 
400 nm/ s [43]. Attached to putative myosin homologues, 
macromolecules could be translocated to the poles in 
~ 20 seconds. This is well within cell division times, so 
this would be a viable polarization mechanism. However, 
no cytoplasmic motor homologue has been identified in 
prokaryotic cells. Furthermore, almost all myosins travel 
towards the pointed ( "+" ) tip along actin filaments [44] , 
so a single polarized MreB bundle would probably only 
support motor-driven localization to one pole while un- 
polarized bundles would not support selective targeting 
to one pole and not the other. 

Static association of proteins to the side of MreB bun- 
dles would, through treadmilling, lead to a translocation 
speed equal to the rate of advection times the monomer 
spacing, Xtreadao, hi the direction of the pointed end 
of the bundle. In steady state, the advection rate 
Xtread = 1.2 s _1 is independent of the buckling force, 
the bundle thickness, or the cytoplasmic concentration 
and leads to v tre ad = aoXtread = 5.92nm/s, shown by 
the dashed line in Fig. 7c. Applying Eq. 7 to a typical 
configuration with p = 1 /itm and L c = 3 /im, this yields 
(ttrans) ~ 12 minutes, which is plausible compared to 
cell-division timescales and comparable to oriC translo- 
cation times [27]. 

Proteins that could bind either directly, or via putative 
tip-binding proteins (analogous to formin in eukaryotes 
[16, 45]) to the barbed end of free protofilaments could 
be translocated in the opposite direction to the tread- 
milling advection. In an unbuckled (Fb — 0) filament, 
lateral protofilaments treadmill at the same rate of the 
net backwards advection of the bundle, accomplishing 
no net movement relative to the cell's axis. In a buckled 
filament, however, the monomer concentration is consid- 
erably above the critical concentration of a free protofila- 
ment, as illustrated in Fig. 3. As a result, unconstrained 
laterally associated barbed ends grow faster than the bulk 
of the bundle, with 

V b = a [Oon C - k otf) - ^tread] 

= a k^, n (c- c c ). (8) 



For each bundle thickness and pitch simulated, the 
translocation times are shown in Fig. 7a and speeds in 
Fig. 7c. Typical bundles with 20000-40000 molecules of 
MreB and 8 — 15 protofilaments thick would transport 
passengers in 4 — 6 minutes, considerably faster than 
laterally-associated proteins and in the opposite direc- 
tion. 

Proteins associated with the slow-growing pointed end 
have a net velocity given by 

V p = O [(K n C - k P oS ) + Xtread] 

= a kP n (c-c c ). (9) 

Note that free pointed ends are disassembling on aver- 
age, though not as rapidly as the treadmilling, so that 
< v p < vtread- These translocation times are shown in 
Fig. 7b and speeds in Fig. 7d. According to Eqs. 8 and 9, 
they are slower by a factor of fcj] n /^ n w 10 compared to 
free barbed ends, taking several hours. These times are 
probably too slow to be biologically relevant. 

These protofilament-associated translocation modes 
offer a non-motor based mechanism for specific target- 
ing of proteins to either pole in cells with a single po- 
larized bundle of MreB. The cell could specify the spe- 
cific pole destination of a particular protein by speci- 
fying which part of an MreB protofilament it binds to: 
barbed-associated proteins would end up at the pole at 
the barbed-tip of the MreB bundle with a speed of v D , 
while laterally-associated proteins would end up at the 
pointed-tip pole with speed v tre ad = o-oXtread- Of course, 
these translocation mechanisms may also supplement a 
(hitherto undiscovered) motor-based mechanism to pro- 
vide targeting to either pole with polarized MreB bun- 
dles. We do not see any way of specific targeting of pro- 
teins to a given pole if the MreB bundle is not polarized 
or if there are anti-parallel bundles, cither with or with- 
out motor proteins. 

If protofilaments dissociate at a significant rate from 
the main bundles then these translocation times repre- 
sent lower bounds. Additionally, any putative MreB- 
binding proteins could strongly affect the polymeriza- 
tion kinetics. For example, ADF/cofilin in eukaryotes 
increase k^ s for actin by <~ 20 times. A bacterial ho- 
mologue of such a protein would decrease the delivery 
time of pointed- associated proteins by ~ 20. Similarly, 
in the presence of profilin, formin increases the barbed 
growth rate of actin by 10- to 15-fold [46], and such a 
modification would decrease the delivery time of barbed- 
associated proteins to within a minute. 

Proteins associated with pointed or barbed ends of 
protofilaments could be translocated towards those poles 
at v p and v , respectively. However the proteins could 
also be directly recruited to distinct poles of the cell, 
due to the free ends of the MreB bundles. The rela- 
tive magnitude of translocation vs. direct recruitment is 
dependent on the number of barbed or pointed protofila- 
ment ends. Tip-associated translocation requires a signif- 
icant number of laterally associated protofilamcnts tips 
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to maximize the translocation flux, while tip-directed po- 
lar recruitment requires unbroken protofilaments to mini- 
mize non-polar binding sites. Observations in C. crescen- 
tus [10] indicate that protofilaments are short, support- 
ing tip-associated translocation as a viable mechanism in 
vivo. 



D. Recovery dynamics of helices 

In experiments on B. subtilis involving fluorescence 
recovery after photobleaching (FRAP) of fluorescence- 
tagged Mbl, four helical turns on one side of the cell's 
longitudinal axis were bleached while the other half con- 
tinued to fluoresce. It took approximately 8 minutes for 
the bleached halves to recover fluorescence at the same 
level as the unbleached halves [15]. We obtain an upper 
bound for this time by calculating the average time un- 
bleached parts of the MreB protofilaments take to tread- 
mill into the bleached regions; one half-turn in the po- 
larized model and one quarter-turn in the non-polarized 
models. Thus, in silico, 

ttR c n^/p 2 + Rl 

tFRAP ~ 7—777 = 7 (10) 

aoSmUAtread O-oMread 

for polarized array structures and half of that for non- 
polarized array structures. For p w 0.5 ^m this yields be- 
tween 3 and 6 minutes. The agreement of this timescale 
with experiment suggest that monomer renewal by ex- 
change with the cytoplasm may not dominate the FRAP 
recovery time. 

Several experiments have applied the MreB-specific 
small molecule A22 to quickly and reversibly break 
down the MreB cytoskeleton by blocking polymerization 
[23, 25]. Cells remain viable after recovery from A22- 
induced disruption of MreB and reform their helical pat- 
terns in less than 1 minute for C. crescentus [25]. In the 
context of our model, recovery from A22-treatment cor- 
responds to re-establishment of the steady-state MreB 
helix from a pool of cytoplasmic monomers. We simu- 
lated our stochastic model from a nucleus of n protofila- 
ments, each of length 3 (as suggested for actin [47]), un- 
der the assumption that the nucleation time is short [31]. 
Thicker bundles (larger n) reach their final steady-state 
length much faster than thin bundles, due to the pres- 
ence of more free filament ends. The equilibration times 
vary between 1 — 5 seconds. Fewer bundled protofila- 
ments took longer to reach a steady state, the longest 
being n = 1, with N = 10000, in which the final length 
is reached in 5 seconds. These are consistent with A22 
recovery timescales. 

We can also address the timescale of breakdown. As- 
suming that A22 simply blocks polymerization, but does 
not change depolymcrization dynamics, then it will take 
t\ = Lfu/(koff + k ff) w 13 minutes for each protofila- 
ment to disassemble with no internal free ends. If each 
protofilament has length L proto with m = Lfu/L proto 
free pointed (or barbed) ends, then we would expect the 



disassembly to be correspondingly faster by t m — t\/m. 
Experiments constrain the actual disassembly time in C. 
crescentus [25] to be < 1 min, implying to > 13 and 
Lproto ~ 300 nm. Analysis of single-molecule experi- 
ments in C. crescentus estimates L pro to ~ 392 ± 23 nm 
[10] by assuming that protofilaments treadmill in place. 
If protofilaments also advect along the bulk MreB cable, 
as we assume in our model, this experimental estimate 
represents an upper bound and is consistent with our re- 
sult. 

These estimates of assembly and disassembly 
timescales of the MreB helix can also be applied 
to the reported midcell condensation of MreB in E. 
coli [12] and C. crescentus [13], and of MreBH in B. 
subtilis [42]. Two possible mechanisms are an elastic 
compression of the intact MreB helix to midcell driven 
by some (posited) motor protein, or disassociation of 
the MreB helix and (transient) association with some 
midcell binding partners. These mechanisms may occur 
in tandem. The timescales of condensation and recovery 
are w 30 minutes in C. crescentus which is considerably 
longer than A22-induced breakdown. The maximum 
filament-end force needed to compress a typical MreB 
helix in our model from L c to a midcell spiral is 
pmax i0p./V, and is comparable to the force generated 
by the RNA polymerase [23] and the anchoring forces 
of integral membrane proteins [48]. Thus, both motor 
driven compression and depolymcrization arc plausible 
mechanisms for the observed midcell condensation. 



IV. DISCUSSION 

Our model explains how polymerization forces in a 
bundle of MreB protofilaments can maintain a helical 
configuration in mechanical equilibrium. Helices are ob- 
served increasingly often within bacteria [4] and our 
model suggests an elastic mechanism for maintaining 
such helicity when filament bundles extend the length 
of the cell. For MreB, our model predicts that the ob- 
served pitches and protein abundance are consistent with 
bundles of 10 — 20 protofilaments thick. 

Helical pitch is directly observable through fluores- 
cence microscopy, and the total abundance of MreB can 
be controlled by an inducible promoter [15]. Increasing 
the MreB expression level should systematically decrease 
the pitch according to the curves in Fig. 2. We have as- 
sumed a cytoplasmic pool of monomers, while there is 
evidence of an oligomeric pool of Mbl in B. subtilis [15]. 
Since only monomers will contribute to force-generation 
through polymerization, an oligomeric pool that was not 
part of the MreB bundle would simply lead to a shift 
along the N axis in Fig. 2 and 3. 

Filament treadmilling is ongoing in the steady-state. 
A significant result is that parts of the filament may 
travel faster than the bulk treadmilling rate, even in the 
opposite direction, while other parts may travel much 
slower. Since the cytoplasmic concentration is held above 
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the critical concentration of free protofilaments by the 
constraining forces at the filament tips, the barbed end 
of protofilaments that associate laterally with the MreB 
bundles can travel at speeds greater than the bulk advec- 
tion speed. These speeds are fast enough to be biologi- 
cally relevant in the transport of polar-targeted proteins 
and the origin-proximal part of the chromosome. These 
transport modes may explain how MreB cables transport 
passengers without motor proteins, and may also rec- 
oncile experimental observations of very different speeds 
of MreB protofilaments: fast speeds of 0.07^iin/s [11] 
may correspond to the lateral protofilaments described 
above, while slow speeds comparable to actin tread- 
milling speeds [10] may correspond to protofilaments in- 
volved in the bulk of the cables. It also suggests that 
observations of bidirectional transport along MreB ca- 
bles [10] does not necessarily rule out a single polarized 
MreB bundle. 

It has been suggested that cell polarity arises from the 
inherent polarity of MreB polymers [16, 22]. The differ- 
ent modes of transport in our model provide a mecha- 
nism for translating the polarity of a cable to the en- 
tire cell. We propose that proteins destined for one 
pole associate laterally with the cable, and are natu- 
rally translocated towards the "pointed" pole. Proteins 
destined for the other pole bind to the barbed ends of 
free protofilaments and are naturally translocated to- 
wards the "barbed" pole. We have demonstrated that 
such translocation happens within a physiologically rel- 
evant timcscalc. Such selective targeting requires polar- 
ized MreB filaments. Proper polarization of the MreB 
cables themselves may require an upstream nucleator, for 
which TipN in C. crescentus is a candidate [49]. While 
antiparallel MreB cables can support transport to both 
poles, we do not see how it can support targeted trans- 
port to specific poles. 

Evidence for polarity of MreB in C. crescentus is mixed 
[10, 22], however both PleC and DivJ appear to have 
targeted localization to specific poles during the cell cycle 
[23] . This is consistent with a single polarized MreB cable 
in C. crescentus. We are not aware of evidence about 
MreB cable polarization in E. coli. While FRAP recovery 
of Mbl cables in B. subtilis is symmetric [15], this could 
be due to lateral exchange and/or the motion of laterally 
associated protofilaments. 

Any transport along treadmilling MreB bundles would 
require either direct association of targeted proteins with 
MreB, or association through an intermediary. Interac- 
tions of only a few cytoplasmic proteins with MreB have 
been identified. These include SetB, implicated in chro- 
mosome segregation [24], RNA polymerase [23], and the 
chaperonin GroEL [23]. Of these, GroEL is known to in- 
teract with a large number of cytoplasmic proteins [50], 
while RNA polymerase is needed for their transcription. 
GroEL is localized to polar and septal regions in E. coli 
[51] while RNAP is associated with nucleoids in punctate 
patterns similar to those seen for helically distributed 
proteins [52, 53]. It would be interesting to see if any of 



these MreB associated cytoplasmic proteins exhibit pole 
directed transport within bacterial cells. 

Transport speeds may be decreased if short protofil- 
aments necessitate multiple binding and release events 
to translocate from one pole to the other. Tip-directed 
transport also allows for direct recruitment to free 
protofilament tips, such as at the bundle ends. The de- 
tails of bundle ultrastructure and of protein translocation 
will be required to sort out these competing effects. 

We have shown that many of the dynamical phe- 
nomenon associated with MreB, including polar protein 
localization, MreB helix dissolution and reformation, and 
FRAP studies, could result from polymerization and 
treadmilling of bundled MreB protofilaments without 
motor proteins. In order to make precise quantitative 
predictions, direct measurements of the kinetic rate con- 
stants for MreB polymerization would be invaluable. The 
Thermotoga-scaled MreB parameters that we used are 
consistent with experiment and require the least amount 
of ad hoc parameter adjustment. Pitch-abundance and 
concentration-abundance relationships did not change in 
character between eukaryotic actin and our Thermotoga- 
scaled MreB kinetics. Conversely, translocation times 
and timescales for dissolution and reformation of MreB 
helices are strongly parameter-dependent. However, our 
qualitative result of bidirectional transport is general 
since, from Eq. 8, is always in the opposite direction 
from the bulk treadmilling. 

There are two significant sources of stochastic effects in 
our model. Stochastic fluctuations due to the number of 
force-bearing protofilaments, n t i P , significantly increases 
the net force applied to the bundle (Fig. 4A). There are 
also large fluctuations in the cytoplasmic monomer con- 
centration in the steady-state, c (Fig. 3), and these sys- 
tematically decrease the net forces (Fig. 4B). We have 
shown that the stochastic effects of c, due to the small 
cell volume combined with the low critical concentration, 
and the stochastic effects of n t i P , due to the relatively 
small number of protofilaments, are both significant. 

We have assumed that MreB polymers are naturally 
straight, without intrinsic curvature or twist. There are 
suggestions that MreB helices are always right-handed 
[7]. Handedness could result from nucleation conditions 
or constraints at the helix ends, which would not affect 
our results. Handedness could also arise from a non- 
zero intrinsic twist r Q in Eq. 4. For small |r | <C i2~ , 
the elastic force in Eq. 5 would remain unaffected. It 
has been hypothesized that MreB exhibits an anisotropic 
membrane affinity requiring one face of the monomers to 
always orient towards the membrane together with signif- 
icant intrinsic curvatures [54]. If so, our polymerization 
dynamics could still be applied though with a more com- 
plex energy functional [54]. However, we would expect 
qualitatively similar results if helical extension remains 
force- limited. 

We have demonstrated that many of the dynamical 
phenomena associated with MreB including polar pro- 
tein localization, MreB helix dissolution and reformation, 
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and FRAP studies, can be addressed by polymerization 
and treadmilling of bundled MrcB protofilaments with- 
out the need to invoke motor proteins. The study of the 
MreB helix ultrastructure, as well as the detailed study of 
the translocation of MreB-associated pole-directed pro- 
teins in the various bacterial species with various MreB 
analogues, will shed further light on the viability of our 
results in vivo. 
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APPENDIX A: GLOBAL ELASTICITY MODEL 

Following Antman [40] , an elastic filament is described 
by the position of its centerline f(s) and an orthonor- 
mal basis of directors {d\(s), (^(s), ds(s)} specifying the 
orientation of its cross-section. This approach has been 
used extensively for DNA [39], but is seldom used for 
eukaryotic actin (though see [55]) since actin does not 
systematically form structures in vivo that are much 
smaller than its persistence length. In contrast, MreB 
forms a helix with a radius of w 400 nm, which is smaller 
than the persistence length of a single protofilament 
(£ p = B/k B T ss 15 /mi) and much smaller than t] p for 
a bundle of several protofilaments. 

The stretching modulus of a filament bundle of actin 
is nE :=y 40n pN where n is the number of protofilaments 
in a typical cross-section of the bundle. For a typical 
filament with n = 15, the energy scale to stretch the 
filament by monomer addition is AU s t re tch = E ncio ~ 
2 x 10 3 pNnm, whereas a comparable energy scale for 
bending the helix is Ai7 bend = Ba /4:R 2 n 2 w 70pNnm. 
We therefore assume that the MreB filament bundle is 
inextensible and parametrize the filament by its arclcngth 
s € [0,Lfu]. The angular strains in the filament are then 
simply d s di = u x di, u = K\di + K2^2 + Td^ where k 2 = 
the curvature and r is the local twist. We set 
d 3 = d s r, following the standard shear-free assumption 
of biopolymers. 

The centerline of a helix with pitch angle 9 (measured 
from the cell's axis, z) and radius R c is 



, sin \ „ . / sin ( 

r(s) = it c cos ( s ) x + R c sm [ —p^—s | y 



Rc 



R c 



+ cos 6s z. (Al) 



Using the tangent, normal and binormal unit vectors, 
{t(s),n(s),b(s)}, and the Frenet-Serret Theorem [40] we 
obtain the curvature of the helix 



sin 2 6 



K = 



Rc 



(A2) 



and the torsion of the filament's centerline, 

t c = -^-cos#sin#. 
R c 



(A3) 



The directors {di(s), ^(s)} are, in general, a rotation of 
{n(s),b(s)} through an angle <fi(s). The difference be- 
tween the physically relevant total twist r and the cen- 
terline torsion t c is tl = t — r c = d s cf> [56]. 

Using the linearly elastic Hamiltonian from Eq. 4 and 
L c = L fu cos 9 and assuming that tl is free to rotate to 
eliminate r(s), we find 



9H 2B . 2 2 

_ = -_sin0 COS 0<n>, 

dH BLfu . 4 2 

dRc = __ RT Sm ^ } ' 
dH B 

dLj~i ~ ~2Rl 



sm 2 6 (l + 3cos 2 6>)(n 2 ). 



(A4) 
(A5) 
(A6) 



Regardless of any mechanisms holding the MrcB bun- 
dle in a helical configuration, one additional monomer 
must provide an energy of (—aodH/dLfu) to polymerize 
itself to the tip of the longest protofilament (s). This pro- 
vides an estimate for the force acting upon the tip of the 
filament bundle, 



F B = 



j7 B sin 2 #(l + 3cos 2 0)(n 2 ) L fa > L c 
[0 Lfu < L c , 



where f B = B/2R 2 C w O.OSlpN. This implies that F B « 
lOpN for n w 10. This is comparable to pull-out forces 
of integral membrane proteins [48]. 

Eq. A4 provides an estimate for the force generated 
by the MreB cytoskeleton against the cell's end caps. An 
estimate for the difference in longitudinal spring constant 
of the cell with and without a properly formed MreB helix 



is, 



Afc 



d 2 n 



2B 



cell 



dL 2 LJ& 



cos6»(3cos6i- l)n 2 . (A8) 



For p = 
For E. 



1 /mi and n = 15, this yields 4 x 10~ 4 pN/ nm. 
coli, the spring constant of the entire cell is 
w 10 3 pN/ nm [57] and is much larger for B. subtilis, thus 
the spring constant differential provided by the MreB he- 
lix is insignificant. Similarly, Eq. A5 provides an estimate 
for the radial line pressure the helix exerts on the lateral 
walls: 



Ff B . 4 2 

7 — ~ ^3 sm 9n ■ 
Lfu i? 3 



(A9) 



For p = 1 /mi and n = 15, this yields 0.03pN/nm. 
Eq. A9 represents a force contrast between pushing the 
cell wall directly above the helix and elsewhere. A bun- 
dle with n — 15 has approximate thickness aQ^/n w 
20 nm and the local excess pressure over the cables is 
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~ 5 x 10~ 3 pN/ nm 2 . In comparison, turgor pressure in- 
side E. coli has been measured at sa 0.1pN/nm 2 [57], 
indicating that the rigidity of the MreB bundles docs not 
directly provide significant structural support for the cell 
well. 

We may also use the elastic model to check the 
self-consistency of our quasi-static helical configuration, 
Eq. Al. At physiological temperatures, all cellular com- 
ponents undergo thermal fluctuations of the order fe^T. 
However, a typical helical bundle with n = 5 and p = 
1 um has an elastic energy of U — 3.5 x 10 5 pNnm ^> 
fcgT = 4.1pNnm. We therefore assume that the heli- 
cal pitch does not significantly fluctuate due to thermal 
effects. An estimate for the timescale for elastic reorga- 
nization of the MreB helix can also be obtained from the 
relaxation time of the lowest hydrodynamic mode in the 
filament [58], 



^elastic 



2 6 Vj4u 
(37r) 4 ln2 B 



0.02 s, 



(A10) 



where rj is the cytoplasmic viscosity. In comparison, the 
characteristic time for polymer elongation is given by the 
steady-state treadmilling rate, 



tpoly 



1 



on 1 on 



A 



tread 



UP Lb _ Lb LP 
^off^on %ff ft ' on 



1.7s. 



(All) 



Since t e i as uc *C t po i y we can assume that elastic relax- 
ation is fast compared to the polymerization dynamics 
of interest. 



APPENDIX B: MEAN-FIELD STEADY-STATE 

By neglecting stochastic fluctuations, we can analyti- 
cally relate pitch, total monomer number and cytoplas- 
mic concentration. The two essential ingredients are sta- 
tionary filament lengths due to treadmilling in the pres- 
ence of elastic constraint forces, and conservation of the 
monomer number. If n t i p G [1, n] protofilaments reach 
each tip, then from Eqs. 1, 2, and 5, 

F B = fo In (j^J n tlp = f B sin 2 6 (l + 3cos 2 #)(n 2 ). (Bl) 

We also have conservation of total number of monomers, 

L r 1 



N = N A c + (n) 



do COS0' 



(B2) 



where N a is the total number of monomers, N A is the 
number of molecules per fiM for a given cell volume, and 
(n) is averaged along the length of the helix. Together 
these equations give a relationship between concentration 
and pitch angle, 



In I ■ | n U p = ^-sm 2 6 (1 + 3cos 2 £ 



fo 



/ | (N - N A cf cos 2 e ) . (B3) 



and an equivalent relationship between bulk thickness 
and pitch angle, 

lli \N A c c {n) a N A c c cos 8 ) nup 

= ^sm 2 d (1 + 3cos 2 0)(n 2 ). (B4) 
Jo 

The mean-field force-concentration relation obtained 
from Eq. Bl and Eq. B4 is plotted in Fig. 4 for n t i P = n, 
which maximizes Fb for the bundle. 



APPENDIX C: STOCHASTIC BUNDLE TIP, 

Nbulk = 2 

If a force F is applied to both ends of a protofilament, 
both k^ n and fcg n will be reduced by a factor of e = e~ F /f° 
[30, 34]. Though the treadmilling rate is unaffected by 
the forcing, the treadmilling concentration becomes 



c = 



off 



Lb 
"-on 



,F/fo 



(CI) 



In the absence of forces, every protofilament of a bun- 
dle treadmills with equivalent rates and critical con- 
centrations. However, force generation by a bundle in 
which the load is applied only on the most advanced 
ntip protofilaments behaves differently Lagging filament 
tips grow without force retardation and with a reduced 
treadmilling concentration, thereby catching up to the tip 
more often than otherwise — increasing n tip and reduc- 
ing the load per loaded filament. This leads to a larger 
total load for the bundle for the same monomer concen- 
tration c. A similar effective increase in force generation 
due to stochastic fluctuations has been noted in passing 
[59], although not analyzed in detail. Here we explicitly 
work out the details for n = 2. For n > 2, the effects of 
force generation by a stochastic n t i P are expected to be 
even more significant, as is seen in Fig. 4. 

We parametrize the system of two filaments by i € 
{0, 1, ...}, the number of monomer spacings between the 
two tips. For any i, two competing Poisson events could 
increase the spacing: addition at the leading tip and dis- 
sociation at the lagging tip. Reduction of i occurs by the 
complementary two events. The master equations lead 
to the following equations for stationary probabilities Pi\ 

i>2: Pi-i(k oS + k on ce) + p i+ i(k on c + k oS ) 

= Pi(k tt + k on ce + k on c + k oB ), 



i = 1 : 2p a (k off + k on cy/e) + p2(k on c + k s) 

= pi(2fc off + fc on c(l + e)), 

i = : 2p a (k oS + k on cy/e) = pi(k on c + k s). 
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The factors of y/e arise when i — and there are two 
leading tips sharing the load, halving the force per fila- 
ment. These equations lead to recurrence relations for pi. 
Although i cannot exceed the length of the filament, the 
Pi vanish exponentially as i increases so we approximate 
imax — oo. Imposing normalization, J^Pi = 1> we nn( l 



P b o = 



C(l - 6) 



(C2) 



where c b = k b ofi /k b n . An equivalent expression for p% 
applies for the pointed end. The average polymerization 
velocity for these stationary probabilities is 



(2k b on cV~e) P v + (k b on 
(k b a ce - • 



ce-k b s )(l- Po ) (C3) 
K s ) + (k b oa c(2V~e-e) + k b oS )p Q , 



and is a nontrivial function of both c and F (through 
e). To solve for the treadmilling concentration at a given 
force, we insert the barbed and pointed equivalent ver- 
sions of Eq. C3 into A p = — A fc , which is a cubic polyno- 
mial in c, and we extract the one real, stable root. 



The relation between the total applied force, Fb, and 
the cytoplasmic monomer concentration c is shown by 
the solid line in Fig. 4. Note that this derivation (and 
the simulations in Fig. 4) imposes a static c, though nu p 
fluctuates. The disagreement with the mean-field result 
(dashed line), indicates that the nu p fluctuations are im- 
portant for the total force. 
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FIG. 1: Five models for the MreB ultrastructure, each com- 
posed of polarized protofilaments. A) Single polarized bundle 
of protofilaments with lateral interactions that prevent rel- 
ative slipping within the bundle. As indicated by the star, 
data in subsequent figures are for this model. B) Polarized 
bundle of protofilaments that freely slip relative to each other. 
C) Unpolarized bundle of protofilaments, with no lateral slip- 
ping. D) Two antiparallel bundles of polarized protofilaments 
that slip relative to each other, though protofilaments within 
a given bundle do not slip. E) Unpolarized bundle of protofil- 
aments that all freely slip with respect to each other. 
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FIG. 2: Steady-state pitch p vs. total molecule number No for 
various filament bulk thicknesses n as predicted by stochastic 
simulation. The mean-field plot is indistinguishable from this 
due to the low cytoplasmic concentration. The rectangles 
represent approximate regions of experimental relevance from 
[7, 9, 14]. In B. subtilis, if the three isoforms bundle together 
into a non-slipping triplex structure, the number of monomers 
will be the sum of each homologue. 
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FIG. 3: Steady-state cytoplasmic concentration c of MreB 
monomers vs. total monomer number iVo from stochastic 
simulation for n protofilaments, showing n = 1, 5, 10, 15, 
20, and 25. The dotted line represents the critical concentra- 
tion. The dashed line illustrates the large standard deviation 
of stochastic fluctuations in steady-state, for n = 5. For all n, 
relative fluctuations were 50% ± 5% and absolute fluctuations 
were within 0.2nM of those shown for n = 5. 
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FIG. 4: Maximal steady-state force generation F by a fila- 
ment bundle with a stochastic a function of steady- 
state concentration c. Various average bundle thicknesses n 
are shown. The dashed lines represent the mean-field pre- 
dicted force-concentration relation (Eq. Bl) if nu v = n inde- 
pendent tips were all sharing the load, (a) We hold c fixed and 
allow an arbitrary bundle stiffness. The solid line represents 
the analytic prediction for n = 2 with a stochastic n t i V . The 
points indicate stochastic simulations, only allowing nu P to 
vary. Fluctuations in nu v allow a significantly increased force 
compared to the mean-field results, (b) The points indicate 
fully stochastic simulations, for the same n as in a), and the 
average c is used. A specific bundle elasticity is imposed by 
forcing F = Fb- Fluctuations in c systematically decrease 
the bundle force compared to a), and this effect is stronger 
for smaller c. 
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FIG. 5: Pitch p as the cell length L c grows for various filament 
bulk thicknesses ribulh as predicted by stochastic simulation, 
for (a) iVo = 3 x 10 3 L C so that an average 3/xm cell contains 
9000 MreB monomers, corresponding to regimes towards the 
low-iVo end of Fig. 2, and (b) No = 10 4 L C , corresponding to 
the high-TVo end. For most bundle thicknesses, the pitch is 
effectively constant as the cell elongates. However, for thicker 
bundles in the low-iVo regime, the helical pitch exhibits a 
significant dependence on cell length, expanding as the cell 
doubles in size. 



FIG. 6: Schematic of the different possible modes of transport 
along MreB bundles. Velocities are absolute with respect to 
the cell in the indicated directions. Bulk treadmilling advects 
bundles that span the cell length at a speed v trea d towards 
the slow-growing pointed ("— ") end. Side associated protofil- 
aments, without constraint forces at their tips, have differ- 
ent polymerization rates and so are not simply advected with 
the bundle. Their tip velocities are v p and vt for pointed 
and barbed ends, respectively. As discussed in the text, vt 
is always opposite the treadmilling direction. Putative mo- 
tor proteins would probably be polarized and would have a 
characteristic speed v mo t- The figure illustrates one polarized 
protofilament bundle, however it is possible that two oppo- 
sitely polarized bundles would exist within a cell — in which 
case the polarity and velocities of the second bundle should 
be opposite the first. 
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FIG. 7: Mean translocation times (ttrans) for macromolecule 
passengers being transported by MreB side protofilaments 
treadmilling on the side of the main cables vs. No. Con- 
straint forces at the ends of the helix keep c above c c , so that 
significant axial movement can be seen with respect to bulk 
treadmilling. Different modes of transport along the cables 
are depicted schematically in Fig. 6. Barbed ends move an- 
tiparallel with respect to bulk advection (a,c). Pointed ends 
move in the same direction as bulk advection, however the 
slower polymerization rates lead to much longer translocation 
times (b,d). Lateral traffic associated with parts of the bulk 
cables move towards the pointed end at ao\ tr ead, indicated 
in c) by the dashed line. Any passenger dissociation and re- 
association from protofilament tips will increase translocation 
times and decrease effective speeds. 
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Symbol Value 


Meaning 


Reference 


Eukaryotic actin kinetic rates for polymerization 




k b 

"-on 


12/iM _1 s _1 


Barbed-end addition rate constant 


[29] 


k p 
™on 


1.3 fjM^s' 1 


Pointed-end addition rate constant 


[29] 


K oS 


1.4s- 1 


Barbed-end dissociation rate 


[29] 


fc off 


0.8s- 1 


Pointed-end dissociation rate 


[29] 


Cc 


0.167 fjM 


critical concentration 


(Eq. 2 and [29]) 


Atread 


0.58 s" 1 


free treadmilling rate 


(Eq. 3 and [29]) 


Assumed MreB kinetic rates for polymerization 




k b 

"'on 


1360 fiM' 1 s" 1 


Barbed-end addition rate constant 


(scaled to fit c c and A) 


^on 


150 fiM^s' 1 


Pointed-end addition rate constant 


(scaled to fit c c and A) 


k b oS 


2.9s- 1 


Barbed-end dissociation rate 


(scaled to fit c c and A) 


h v 
K oS 


1.7s- 1 


Pointed-end dissociation rate 


(scaled to fit c c and A) 


C c 


0.003 fjM 


critical concentration 


[31] 


Atread 


1.2s- 1 


free treadmilling rate 


[10] 


B 


1.0 x 10 4 pNnm 2 


Bending modulus of MreB protofilament 


(from actin, [60]) 


E 


20 pN 


Stretching modulus of MreB protofilament 


(from actin, [60]) 


ao 


5.1 nm 


Monomer spacing of MreB 


[31] 


Rc 


400 nm 


Radius of B. subtilis cell 


[14] 


L c 


(2 - 4) x 10 3 nm 


Length of cylindrical portion of B. subtilis cell [14] 


V 


10" 9 pNs/nm 2 


Cytoplasmic viscosity 


[61] 



TABLE I: The parameters used in this paper for MreB and its paralogues. Only c c and ao have been directly measured for MreB 
from T. maritima. Other assumed parameters have been extracted from the properties of eukaryotic actin, which is homologous. 
The kinetic rates chosen for MreB recover c c [31] and the in vivo treadmilling rate [10], but leave the pointed/barbed asymmetry 
unchanged. The elastic moduli are for a single protofilament of F-actin (n — 1). 



